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ABSTRACT 

The Aeronautics Department shock tube has been developed and 
instrumented. The shock tube has been used in an experiment in which 
the thermal conductivity of argon has been determined in the temperature 
range of 1500 - 5000 °K and at relatively high pressures, 10 to 30 
atmospheres . 
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CHAPTER I 



INTRODUCTION 

Although the shock tube has been known as an experimental device 
in the study of gas dynamics for the past 70 years^ since Paul 
Vieilli in France (1899) made first known use of it, it was not until 
the early 1940' s when better methods and equipment were developed to 
measure the phenomena occuring in the gases tested, that it began to 
enjoy much application or widespread use. 

With the rapid technological development of electronic measuring 

equipment permitting the use of fast response high sensitive instruments 

the shock tube has become a very useful experimental tool for study 

and research in the fields of gas dynamics. Specific applications include 

the study of shock waves, shock fronts and interactions, aerodynamics 

for the study of the subsonic, transonic and supersonic flow regions, 

as well as a means of driving hypersonic wind tunnels^, chemistry 

F21 

for the study of high temperature gas reactions 1 - , physics for study 
of conversion of energy from one form to another, i.e., translational 

m 

to energy of internal vibration on rotation of molecules 1- , in the 
field of planetary-entry fluid physics connected with the design of 
the spacecraft's heat shield^ and many other applications. 

The purpose of the work presented in this thesis, was first to 
install instrumentation and leave in working condition a newly 
developed shock tube for the Department of Aeronautics of the Naval 
Postgraduate School and after that was accomplished, use it to 
determine the thermal conductivity of a noble gas (argon in this 
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case) in the range from 1500 to 5000 degrees Kelvin at relatively 
high pressures i.e. 10 to 50 atmospheres. Although several authors 
have reported results in this type of experiment^ , they all 
have been done at low pressure i.e. around 1 atmosphere and data for 
the same experiment worked here was not available. The value obtained 
here, for the calibration parameter of the shock tube, g , is some 
15% higher than similar parametric values verified in previous shock 
tube experiments^ , but the initial conditions to determine it, 
where again different, in this case, the calibration runs were 
accomplished in the range from 1 to 10 atmospheres of pressure. 

As a preview of the following work. Chapter II describes the 
equipment, instrumentation and procedures used in the course of this 
work. Chapter III is related to the general theory behind shock tube 
operation as well to the theory related to the experiment accomplished. 
Chapter IV is devoted to the experiment itself, the data reduction 
and presentation of results, and Chapter V gives the conclusions of 
the present work. 

It is proper to hereby acknowledge the kind assistance of several 
people who contributed to the completion of this work. The author 
wishes to express his appreciation to his advisor. Dr. Daniel J. Collins, 
for his guidance and assistance during the entire time this work was 
carried out; to Mr. Dana Maberry and Mr. Norman Leckenby for their 
cooperation and assistance in the performance of the experimental 
program; to Mr. Robert Besel for his assistance with the graphics work. 

Last, but not least, the author is particularly indebted to his 
wife, Maria Angelica, for her assistance in the preparation of the 
manuscript and for her understanding and patience during the entire 
course of studies. 
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CHAPTER II 



DESCRIPTION OF EXPERIMENTAL EQUIPMENT 
I. SHOCK TUBE 



General 

The shock tube shown in Plate I was constructed at the Naval 
Postgraduate School, Monterey, California. The equipment has been 
installed in one of the Aeronautical Laboratories of the School. 

Shock tube configuration was chosen with two definite purposes 
in mind, namely to be used as a laboratory equipment for the courses 
in the Gas Dynamics option in the Aeronautical Engineering curricula 
at the School and at the same time as a research aid in the field 
of high temperature gas dynamics. 

Essentially a shock tube consists of a rigid cylinder (although 
different cross-sectional shapes are possible), divided into two 
sections by a thin diaphragm mounted normal to the axis of the tube. 
The principle of operation calls for a pressure difference to exist 
across the diaphragm. The high-pressure side is called the Driver 
section and a low pressure region is called the Driven section. 

Figure 1 shows a schematic of the shock tube sections concerned. 
The material used is a stainless steel tube, 3 in. in internal 
diameter, with a wall thickness of 0.5 inc. and a highly noli shed 
internal surface. The overall length of the tube is approximately 
25 ft. constructed in four sections of different lengths. This 
arrangement facilitated the construction of the tube and at the same 



13 




PLATE I 

SHOCK TUBE AT NAVAL POSTGRADUATE SCHOOL 

MONTEREY, CALIFORNIA 
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SCHEMATIC OF SHOCK TUBE SECTIONS 



time provides flexibility in the case that different driver or driven 
tube lengths are required. The first section of tubing, with a length 
of 7 ft. 3 in. is the driver section. The other three pieces of the 
tube make up the driven section. Between the driver section and the 
driven section is inserted a transition section which holds the valve 
body and diaphragm section. Connecting the three tube sections that 
make up the driven end are two 6 in. outside diameter instrumentation 
sections of 4.5 in. length. These have two instrumentation ports 
at each side of the tube and two on the upper part, intended to house 
the different types of detectors which can be used, such as pressure- 
transducer, shock detectors, temperature-sensitive resistance (thin 
metal film), ionization detectors, light screen schlieren, etc. 

The end of the shock tube has another instrumentation station, 
with four ports located 90 ° apart of each other. The end plate of 
the tube, also has provision for the insertion of instrumentation. 

All instruments and plugs used in the instrumentation sections are 
flush mounted with the inside wall. 

Different sections are bolted together and end flanges and joints, 
are sealed by means of Vi ton "0" rings. The ratio between the length 
of the tube and the inside diameter is 76 and the ratio between the 
length of the driver and driven section is 0.368. 

As indicated before, the driver and driven section are separated 
by a diaphragm. In this shock tube, there are two diaphragms separated 
by a distance of 2.5 in. called the diaphragm chamber, which will 
be described later, permits a very ingenious but simple way of ob- 
taining a rather precisely controlled bursting pressure. The shock 
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tube assembly is supported as shown in Plate I. 



Driver End 

A very important consideration in the operation of the shock 
tube is the vacuum and gas filling system required. Figure 2 shows 
a schematic of the system used for the driver end section of this tube. 
The vacuum pump evacuates the high pressure region of any remaining 
gas, air and moisture before the section is filled with the gas to be 
used. Gas is admitted to the front end section of the shock tube from 
a high pressure bottle of gas by means of tubing through a regulator 
valve connected directly to the bottle and a pressure control valve 
located beside the tube. By means of a "T" connection two high pressure 
bottles of either the same or different gases, can be used in the 
charging of the shock tube. The pressure of the gas in the tube is 
obtained from a pressure gauge located in the upper part of the front 
end section of the tube. To provide a means of evacuating the gas, 
after the shock tube has been fired, a bleed valve is located in the 
bottom part of the front end section. Plate 2 shows the actual set-UD 
of the driver end section. 

Driven End 

The vacuum system of low pressure region or driven end is an impor- 
tant part of shock tube operation. Figure 3 shows a schematic drawing 
of the different parts which make up the vacuum system. Its details 
and purposes are described below. Basically the main concern faced in 
the design of the vacuum system, is the range of vacuum required. 

With this criteria in mind, the selection of a system which will be 
adequate and expeditious, is of prime importance for optimum utilization 
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FIGURE 2 

SCHEMATIC OF DRIVER END SECTION 
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PLATE 2 

VIEW OF DRIVER END 
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FIGURE 3 



SCHEMATIC OF DRIVEN END 
VACUUM AND FILLED SYSTEM 
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of the shock tube. 



Connected directly to the shock tube valve body is the main 
valve, which controls the operation of the vacuum system. The main 
valve is connected to a General Electric, Model 5KC45PG4X, vacuum 
pump which is designed to initially evacuate the driven end and is 

_3 

capable of pulling up a vacuum of approximately 10 mm. of mercury. 

Also connected to the main valve is a diffusion pump, which is capable 
of pulling up vacuums of 10"^ torr.- The diffusion pump used, is 
manufactured by National Research Corporation (NRC), a subsidiary of 
Norton Company of Newton, Massachusetts. The system also consists of 
a Circular Chevron Cryo-Baffle and a pneumatically operated High Vacuum 
Slide Valve, both manufactured by NRC. The equipment is shown schemat- 
ically in figure 3. Two thermocouple vacuum gauges and one ionization 
gauge tube, are used to measure the pressure. The thermocouples 
gauges are manufactured by NRC and the type used is 0531. The vacuum 
ionization control gauge instrument is also manufactured by NRC and 
its model is NRC 831. The ionization gauge tube is a RG 75-K manu- 
factured by Vacuum Electronics Corporation (Veeco), Long Island, New 
York. The vacuum gauge instrument, vacuum ionization gauge instrument 
and the control switches for the diffusion pump, vacuum pump and pneu- 
matic actuator for the high vacuum slide valve are mounted in an instru- 
ment panel for easier control and operation of the vacuum system. 

Plate 3 shows the actual set-up of the driven end vacuum system and 
instrument panel respectively. 
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PLATE 3 

DRIVEN END VACUUM SYSTEM 
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II. INSTRUMENTATION 



General 

The successful use of the shock tube depends entirely on suit- 
able instrumentation for the required flow measurements. Although the 
type of instruments to be used depend largely on the experiment being 
run, accuracy needed, and money available, they all require a precise 
and fast response, since the duration of the shock wave in the tube 
is of the order of milliseconds. Consideration is given here to the 
actual instrumentation used in the shock tube at the Naval Postgraduate 
School . 

Temperature sensitive resistance detectors 

One of the fundamental measurements in shock tube tests is that 
of the shock speed. One way to accomplish this measurement is by 
using temperature sensitive resistance detectors or most commonly 
known as thin film gauges. These devices are of local manufacture, 
using a short length of glass or pyrex about 0.25 in. in diameter, with 
one flat end and edges carefully polished. A film of Liquid Bright 
Platinum No. 05, manufactured by Engelhard Inc., East Newark, New 
Jersey was used and a thin strip was painted over the flat end and 
extended through the sides of the tube; the film should be around 1/32" 
in width and its length and thickness such that the resistance is about 
100 ohms; after that is accomplished it is let to air dry for a minimum 
of one hour, then baked in an oven to 1240°F for 20 minutes and then 
air cooled; the process is completed by annealing three times to 1100°F 
for 1 hour, allowing complete air cooling every time. Leads are 
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soldered to the sides of the gauge by using Eccobond Solder V-91 , a 
product of Emerson & Cuming Inc., Canton, Mass., and after that, the 
gauge is placed on a phenolic plug, for insulation purposes, to hold 
it in place in the instrumentation port plug of the tube, as shown in 
Figure 4. Although "0" rings are provided for sealing, in order to 
avoid any possibility of leak through the thin film and the phenolic 
plug, the free space between them was filled up with Eooxy made with 
APCO 210 resin and APXO 180 hardneer, products of Applied Plastic Co., 
Inc. (APCO), El Segundo, California. 

The response of a thin film gauge operates, dependent on a change 
in resistance with a change in temperature of the thin metal film. 

In the case of the shock tube, this rise in temperature is provided 
by the shock wave. With a constant current passed through the metal 
film, a rapid increase in voltage is obtained as the shock wave passes 
over the gauge. This voltage properly amplified, is fed into an 
oscilloscope and recorded. 

Three thin film gauges were used for the measurement of the shock 
wave speed. Each one of them was mounted in one of the instrumentation 
station ports, provided for that purpose, as shown in Figure 1. 

Calibration of these gauges was not necessary because the signal 
obtained from them, was used to trigger the oscilloscope and from 
their stepwise signal, time resolution was obtained which allowed 
determination of the shock speed. 

Heat transfer thin film gauges 

This type of gauge is the same as the one described above. A 
different name is used for purpose of identification since its intended 
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use is different. The principle of operation is the same as the one 
described for the temperature-sensitive resistance detectors and the 
manufacture process was done in the same manner. This gauge is located 
on the end plate of the Driven end of the shock tube, mounted flush with 
the interior surface of the plate and held in place by means of a phenolic 
plug of similar construction to the one shown in Figure 4. The gauge 
measures the temperature of the end wall before and after the arrival of 
the shock wave. The signal produced by the gauge is fed into an oscill- 
oscope and is recorded photographically. 

It was necessary to calibrate the heat transfer gauge in order to 
establish the temperature coefficient of resistivity (a). The cali- 
bration was done by placing the heat transfer gauge into a bath of 
ice water and recording the initial gauge temperature and resistance; 
the water was then slowly heated to a temperature approximately of 
90°C, and readings of temperature and resistance of the gauge were 
recorded during the process. The following relationship^’^ was 
used : 



r f = r 0 [i + -(t f - y ] 

where 

Rp = gauge resistance 
Tp = gauge temperature 
R q = cold gauge resistance 
T q - initial temperature 

The temperature coefficient of resistivity (a) was obtained by 
plotting the change in resistance versus the change in temperature. 
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FIGURE 4 

THIN FILM GAUGE ARRANGEMENT 
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The resistance of the thin film gauge was approximately 100 ohms, 
and was found to be a linear function of its temperature. Several runs 
were made for the calibration and one example is given in Figure 5. 

The plot was fitted using the least squares method. The temperature 
coefficient of resistivity was found to be 0.00104 ohm/ohm deg. C. 

Electronics 

In Figure 6 is shown a schematic of the electronic set-up used 
in the shock tube. Gauges Nos. 1, 2, and 3, represent the thin film 
gauges used for the measurement of temperature difference at the end 
wall. The gauges are connected to a bridge as shown in Figure 7. 

Since the resistances R-| , R 2 , R3 and R^ are much larger than Rp, the 
current I through the gauge is maintained steady, thus as the change 
in gauge resistance is attained with the passing of the shock front, 
a stepwise signal in voltage is produced and this signal fed into the 
oscilloscopes. In the arrangement shown in Figure 6, the signal from 
gauge No. 1, goes to an Amplifier where it is amplified about two 
hundred times and from there, is used to trigger both oscilloscopes. 

The output from gauges No. 2 and 3 are connected to a dual beam oscill- 
oscope, where the distance between the two peaks in their traces, gives 
the time delay needed to compute the shock wave speed. The output 
signal from the H.T. gauge is connected to the second single beam 
oscilloscope. The electronic equipment used is a Wide band Differ- 
ential D.C. Amplifier Model 885, manufactured by Astrodata, Inc., 
Anaheim, California; a Type 551 Dual Beam Oscilloscope and a Type 549 
Storage Oscilloscope, both manufactured by Tektronic, Inc,, Portland, 
Oregon. 
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FIGURE 5 

TEMPERATURE COEFFICIENT OF RESISTIVITY (oc) 
FOR HEAT TRANSFER THIN FILM GAUGE 
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SCHEMATIC OF ELECTRONIC ARRANGEMENT 
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SCHEMATIC FOR THIN FILM GAUGES 



III. OPERATING PROCEDURES 



Filling and Vacuum System 

The first step in operating the shock tube, is to obtain a vacuum 
in the Driven end. With reference to Figure 3 and starting with all 
valves closed, the vacuum pump is turned on, and when ready to pull 
down vacuum, control valve No. 2 and the Main valve are open; control 
valve No. 1 connecting the vacuum pump with the diffusion pump, is a 
two-way valve, which in the open position, closes the line to the 
diffusion pump. Working in that condition it is possible to obtain a 

_3 

vacuum of approximately 10 mm. of mercury; readings are obtained 
from either the vacuum gauge and/or thermocouple gauges. If a higher 
vacuum is needed, the diffusion pump is used and one opens the vacuum 
slide valve and close control valves No. 1 and 2. In this condition 
the system is operated until the desired vacuum is obtained. Although 

it was not tested during the present work, with the diffusion pump it 

_g 

is possible to obtain vacuums in the neighborhood of 10 mm.' of mercury. 
When the desired vacuum is obtained, the main valve is closed and then, 
by means of the quick acting filler valve the working gas is pumped 
into the Driven end. The gas is usually stored in high pressure 
bottles . 

The next step is to fill the Driver end with the gas to be used, 
to the desired pressure. A vacuum pump could be used in the driven 
end, in order to assure a clean and moisture free atmosphere before the 
working gas is admitted. With reference to Figure 2, the gas, usually 
stored in high pressure bottles, is fed into the driver section, by 
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means of the pressure control valve and at the same time to the dia- 
phragm chamber. Assuming a pressure P^ is desired in the driver end, 
the diaphragm chamber is filled up to a pressure P^ less than P^. 

The chamber, actually a very small compartment of the tube in length, 
is sealed off by two diaphragms, the first being thick enough to 
withstand the pressure difference P^ - P^ and the second one to with- 
stand the pressure P^ but not P^. When the required pressure has 
been reached in the driver section, the diaohragm chamber is either 
rapidly evacuated or expanded and the two diaphragms are then burst 
almost at the same time. This rather simple method of rupturing the 
diaphragm has the main advantage of giving a precise control pressure 
at which the diaphragms are ruptured and permits firing of the tube 
at the desired moment. This ability to control the firing of the 
tube is highly desired if for any reason something went wrong with 
the electronics and data recording equipment at the last moment. 

To evacuate the shock tube after it has been fired, bleed valves 
are provided along the sections of the tube as indicated in Figures 
2 and 3. 

Data recording 

The successful run of the experiments depends entirely in the 
ability of recording the required data. Due to the short time every 
shot lasts, the low output signals of the thin film gauges and the 
high sensitivity of the equipment to be used, this part of the experi- 
ment required some coordination of the many factors involved. 

The data recorded was the following: driver end pressure, 

driven end pressure, initial temperature in the shock tube, time 
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required for shock waves to travel a fixed distance and voltage rise 

in the end wall at shock arrival. The first three measurements were 

taken by direct reading of the instruments provided for that purnose. 

The time required for the shock wave to travel a fixed distance in 

the shock tube, namely the distance between thin film gauges No. 2 

and 3 ^igure 6] was Gained by means of a dual beam oscilloscope and 

recorded using a Polaroid film. A typical trace of the signal obtained 

is shown in Plate 4. With above data, the speed of the shock wave was 

f2l 

determined and knowing the speed of sound in the gas used L J the shock 
wave Mach number was computed. 

The temperature rise in the end plate is given by the following 
relationship: 



AT 




ae 



The temperature coefficient of resistivity (a) was previously 
determined for the gauge as described before; the current (I Q ) and 
resistance (R ) were measured before the runs. The change in voltage 
was recorded from the oscilloscope hook to the heat transfer thin 
film gauge shown in Figure 6. A typical trace of the signal obtained 
is given in Plate 5. 



33 




PLATE 4 

TRACE OF TEMPERATURE SENSITIVE 
RESISTANCE DETECTORS 
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PLATE 5 

TRACE OF HEAT TRANSFER THIN FILM GAUGES 
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CHAPTER III 



THEORY BEHIND SHOCK TUBE OPERATION 
I . GENERAL THEORY 

The theory of the shock tube is well known and clearly described 
in many textbooks. The same comment can be made in regard to the theory 
of shock waves. It is not the intended purpose of this thesis to go 
deeply into the deviation of the different relations governing the 
actual phenomenon occuring in the tube, but to give a' simple presentation 
on the production of a shock wave and then look in more detail at the 
theory governing the experiment which was performed. 

In Figure 8 is shown the thermodynamic history of the shock tube. 
Figure 8a and 8b corresponds to time t = 0. The Driver end corresponds 
to region (4) and the Driven end to region (1). At t = 0 the temperature 
of both regions is the same. 

At the appropriate time the diaphragm is broken. A shock wave (S) 
is propagated into the driven end and a rarefraction (or expansion) wave ' 
(R) is propagated back into the driver end. This is shown in Figure 8c 
at time t-j > 0. As the shock front travels into region (1), the fluid 
is compressed to a higher pressure P 2 (Figure 8d) and higher temperature 
(Figure 8e) in a nonisentropic process. On the other hand, as the 
rarefraction wave travels into region (4), the fluid is expanded to a 
lower pressure and a lower temperature T^ by a process which is 
essentially isentropic. Thus, although the pressures P 2 and P 3 are equal, 
there are two regions between (R) and (S). These two regions correspond 
to regions (2) and (3) in the graphs of Figure 8. The gas in region (2) 
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FIGURE 8 

THERMODYNAMIC HISTORY OF THE 
SHOCK TUBE 
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has been compressed and heated by the shock wave, however, the gas in 
region (3) has been expanded and cooled. The narrow zone called the 
contact surface or contact discontinuity indicated by (C) in Figure 8c 
represents the line where a temperature and density discontinuity exists, 
although pressures are the same on each side of (C). Physically it 
denotes the position occupied by that gas which was originally at the 
diaphragm and is the point where the experimental gas and the driver gas 
make contact. This "contact surface" travels along the tube behind the 
shock front. The above explanation is the ideal one. In the real case 
some mixing of gases is usually present at the contact surface which 
causes the temperature drop to occur in a finite time period. 

The most common notation for the different regions is the one used 
in the above explanation. Another region, called region (5), is present 
in the shock tube operation where the shock wave is reflected at the end 
wall of the tube with a further rise in temperature and pressure exists. 
This region is shown in Figure 9, which actually represents a physical 
plot (distance vs time) of the phenomenon initially produced in the shock 
tube after the diaphragm has been burst. It can be seen in this plot, 
that the rarefraction (or expansion wave) covers a zone characterized by 
falling temperature, pressure and density below the initial conditions. 
These waves spread with time and the fluid properties changes occur in a 
more gradual process as compared with the shock front and contact 
surface. Due to their spreading action, these waves are called the 
expansion fan, with the initial and terminal wave called the head and 
tail waves respectively. 
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ShocK Tube At Time t, >0 




Distance 



FIGURE 9 

PLOT OF DISTANCE VS. TIME IN SHOCK 
TUBE AFTER DIAPHRAGM IS BROKEN 
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Shock Tube End Wall 



The theoretical analysis for shock tube phenomenon is greatly 
simplified by the assumption of one-dimensional flow, where the flow 
properties are considered to be uniform over any cross-section and to 

vary only with distance along the tube and in general with time. 

< 

The fundamental equations governing a non-stationary one-dimensional 
flow are the continuity, momentum, energy and equation of state and they 
are given by the following relationships respectively: 
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The most common procedure to solve for the relationships applied to 
shock tube operation is using the above equations together with the 
methods of characteristics. A detailed coverage of these derivations is 
given in Reference [1] and [2], and for the method of characteristics 
itself in Reference [9]. 
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Equation (5) relates the initial applied pressure ratio across the 
diaphragm to the shock wave Mach number. It gives a measureing of the 



strength of the shock in a given test gas and shows how for a given 
pressure ratio, the shock strength depends on the specific-heat ratio 
and sound speed in the high pressure gas. 

The strongest shocks can be obtained by using a driver gas having 
a high speed of sound and low specific heat ratio. 

Equations (6) and (7) are useful to determine the state of the gas 
behind the reflected shock. In the form presented they are specially 
suited for the experiment, leading to a full solution of the pressure 
and temperature behind the reflected shock in terms of the incident Mach 
number and initial pressure and temperature, parameters measured directly 
from the experiment. 



The region behind the reflected shock wave (region 5) is idealized 
to consist of a hot semi -infinite gas adjacent to a semi -infinite solid. 
Viscouse dissipation may be neglected and the pressure assumed constant 



conditions, we may write the equations of continuity and energy as 



II. THEORY FOR EXPERIMENTS 



in the gas boundary layer near the end wall of the tube^. Under these 



( 8 ) 




and 



( 9 ) 
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The continuity equation is satisfied by introducing a stream 



function \ p, such that: 



( 10 ) 





a^ 

at 



Changing coordinates from x to \p, the energy equation becomes: 



(ID 

where 



1 ( _36_ x _ _J_ / MU x _30_ x 

a w v at ; " 34; v 0 3if; ; 



( 12 ) 



a 



w 



w 



o c 
M w p 




, K(0) 




With the above transformation, no characteristic time appears in 

2 

the equation, leaving 0 as a function of the expression — *4 — , implying 

V 

that the wall surface temperature, T , is a constant. 

w 

Introducing a similarity parameter 



(13) 



4> 

n = 

(2a w t} 



1/2 



The one dimensional time dependent energy and continuity equations near 
the end wall of the tube reduce to 



(14) 



d / Mil x _de_ X 
dn 0 dn 



+ 




0 



subject to the following boundary conditions: 



(15) 



0 ( 0 ) = 1.0 



0 (oo) 
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The subscripts w and °° refer to numerical values evaluated at the 
end wall condition just after the shock wave arrival and to any quantity 
measured behind the reflected shock wave, respectively. 

Equation (14) can be integrated from n = 0 to n = 00 , after the 
temperature dependence of the thermal conductivity is established. 
Reference [5] and [6] assume the thermal conductivity to vary as a power 
law in temperature 




With that assumption, equations (14) and (15) can be used to 
obtain the following relations: 



(17) 



0 co = -r~ = 

w 







n=0 



By equating the heat flux at x = 0 to the heat flux at the surface 

[51 

of the solid end wall it can be shown 1 - J that 



(18) 



'’w 




T w- T 



w 



s 



where 



(19) 



A 




( k s p s 
^ k w p w 




1/2 



In the above equations, the subscripts s and w refer to properties 
of the gauge and properties of the gas at the wall temperature respectively. 
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Since the quantities 0^ and q w can be experimentally evaluated 
from each run, a series of tests yields an experimental relation between 
them for the gas tested. Equation (16) established that the relation of 
these quantities is dependent on the exponent a. Then, by comparing the 
theoretical and experimental data it is possible to choose a value of 
"a" which best fits the data. The numerical value of "a" was computed 
by minimizing the sum of the squares of the deviation between the 
experimental and theoretical curves. Once it was established, the 
temperature dependence of the thermal conductivity is given by equation 
(16). 
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IV EXPERIMENTAL RESULTS AND DISCUSSION 
I. DESCRIPTION OF EXPERIMENT 



Purpose 

The purpose of the experiment described herein was to look at 
thethermal conductivity of argon at relatively higher pressure, i.e. 

10 to 50 atmospheres and to comnare the results obtained with those 
similar experiments, done at/or around atmospheric pressure. Previous 
data for experiments having the same conditions of the one done for 
this thesis was not available, so a comparison in that sense was not 
possible. 

Calibration Runs 

With reference to the theory and equations presented in the pre- 
vious chapter, it was necessary to calibrate the thin film heat transfer 
gauge to obtain a value of B , in order to compute later the experi- 
mental values of q . This calibration was done in the shock tube, 
w 

by making a series of runs (30 in total), where argon was heated to 
temperatures ranging from 1000 to 1500 degrees Kelvin. In this range 
of temperature, the thermal conductivity of argon has been experi- 
mentally determined. To determine the heat transfer rate q w , a value 
of the parameter "a" is required. The value a = 0.697 taken from 
Collins and Menard^ was used. Equation 18 was then used to obtain a 
value for $ s = 0.062060 cal/cm^ deg. K (sec)^ with a standard 
deviation of + 0.010. This value was found to be around 15% higher 
than others previously determined. Due to this reason, and that the 
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calibration runs were accomplished also at relatively high pressures, 
between 1 and 10 atmospheres, the validity of the use of the parameter 
a = 0.697 is questionable. Unfortunately, due to lack of time it 
was not possible to go into the determination of the optimum value of 
"a" to be used in this case. 

Experimental Runs 

The experimental runs were made in the shock tube previously 

described using the data recording equipment and instruments mentioned 

before. The gas tested was Argon using Helium as a driver. The runs 

were made, as to obtained temperatures between 1500° to 5000° K and 

pressures of 10 to 50 atmospheres in the region behind the reflected 

shock wave (region 5). In the determination of thermal conductivity, 

roi 

the value of k was obtained from the values given in Cook 1 . Thus 
the problem was reduced to find a value of "a" for equation (20) to 
fit the data. 

(20) K (e) = k = e a 

A list of data recorded, as well as the values computed is given in 
Plate 6. 



II. DATA REDUCTION 

For every run made, the data obtained was the driven end pressure 
(p-|), the shock tube temperature (T ), the time resolution for the 
shock wave mach number computation and the voltage rise at the end 
wall thin film heat transfer gauge, used for computations of temperature 
rise in the end wall. With that data, it was then possible to obtain a 
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PLATE 6 

TABULAR SOLUTION OF EXPERIMENT 
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dimensionless plot of 0^ = T^ / T w vs (de/dn) w = for the experiment 
accomplished. After that was done and by using equation (14) chapter 
III and by comparing the theoretical with the experimental values 
obtained and minimizing the sum of their deviations a plot of the 
residuals versus "a" was obtained and from it, the best value of "a" 
to fit the data was selected. 

In appendix I are given sample computations indicating the 
breakdown for equation (14) obtained in order to solve it using the 
IBM 360 computer of the school and for the determination of q w 
for the experiments run. 

Although several partial computer programs were written, 

Appendix II shows the final computer program written that solves 
for the best value of a, having as inputs the initial data obtained. 

The output of those programs was divided in three steps; the first, a 
tabular array given for each run the initial data and the computed 
values obtained for the mach number, temperature raise in the end wall, 
pressure and temperature in the region behind the reflected shock, 
final temperature in the end wall, a value of 0 = T/T and a value for 
q w> With these last two quatnities, a subroutine called W4PL0T was 
used to obtain a graphical representation of the experimental results. 
Then, equation (14) was solved, using a subroutine called SYSZDE and 
from it, a theoretical value of 0 = T /T was obtained, for different 

00 OO Yl * 

values of "a". The residuals between the theory and experiment were 
then obtained and these values, plotted versus a, yielded the selection 
of the optimum value of "a". This was represented in a graph, using 
again the subroutine VJ4PL0T. The final plot obtained is shown in 
Figure 10. 
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Residual 




FIGURE 10 

PLOT OF EXPONENT a VS. RESIDUAL 
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With the computer program shown in appendix III, a graph of the 
experimental data, together with the curve for best fit was then ob- 
tained and at the same time, theoretical curves for other values of 
"a" chosen arbitrarily were plotted for comparison, as shown in Figure 
11 . 

Due to the relatively low time resolution used in the oscillo- 
scope connected to the heat transfer thin film gauge, 1 millisecond 
vs 10 microseconds of reference 6, the flat plateau at the beginning 
of the trace, as the one mentioned in that reference, was not visible 
in all the pictures. 

In approximately one third of the data a flat plateau was iden- 
tifiable. Another reference point visible on all the traces was used. 
This plateau point was some 25% higher than the reference point. A 
correction factor of 1.25 was then used for all the data. 
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CHAPTER IV 



CONCLUSIONS 

The experimental data is preliminary in nature, nevertheless 
the thermal conductivity of Argon at higher pressures 10 to 50 atmospheres 
has been determined. In the relationship given by equation 16, the 
experimental value of a = 0.4 was obtained with k w computed as shown 
in appendix II, step AKW. No previous data in this range of pressures 
has been obtained so that it is not possible to make comparisons. 

It is necessary perhaps, in some future study of the thermal 
conductivity of Argon, to explain from a theoretical view point the 
decrease in the value of "a" with pressure, from that obtained at 
lower pressures. 

It would appear that it is also necessary to repeat some aspects 
of the experiment in order to verify some of the results obtained. 

This applies particularly to the value of 6 $ obtained. 

The purpose of this thesis was to develop a fully functioning 
shock tube and to demonstrate the functioning of the shock tube 
by means of an experimental problem. 

The shock tube is now a working tool in the Aeronautical Engin- 
eering Department and some preliminary results have been obtained 
in an experiment to determine the thermal conductivity of Argon in 
the pressure range from 10 to 50 atmospheres. 
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APPENDIX I 



SAMPLE CALCULATIONS 



a) Solution of equation (14) for data reduction: 



(14) 



K(e) de 

dn l 0 dn 



+ n 



de 

dn 



from equation (12) Chapter III is found: 
K(6) = 0 a 



Replacing in equation (14): 



,a-l d‘e 



dn 



d / n a-l\ d0 
"dn (e } 



+ n 



de 

dn 



= 0 



reducing above expression: 



( 21 ) 



d 2 e , a-1 

dn 2 



e 



de 

dn 






Equation (21) represents a second order nonlinear differential equation. 
In order to solve it, it was reduced to a system of two first order 
differential equations by the following transformation: 



Let 



x-j = e 



de 

dn 



Then, the final equations are: 
(22) x, = x 0 



(23) 



x ] 3 X 2 " nX l ] 3 x 2 
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b) Determination of q for experiment runs: 

w 



The value of B $ was known to be 



B s = 0.06026 



cal . 



k cm^ K°(sec)^ 



From equations (17), (18) and (19) Chapter II, is obtained: 

de 



n w 



dn ‘ ]j 



T e s V T s 



71 B w T w 



let 



AT = T - T 
w s 



Using equation of state for a perfect gas: 



T 3 = T 

w w W' 



c = T 
w p w 



k P c / c 

w w p _ I k p T _e_ 

R T V W w w R 

W 1 



For above expression, the following holds 



P = P in units of [cal /cm" 5 ] 

w 00 

k is given in units of [ — — — ] 
w cm°K sec 

T is given in [°K] 

3 

P = (P mm Hg) (3.184 10' 15 ) is a conversion factor 

c 

= 2.5 for Argon 

K 



Then: 



n w 



d6 

dn 



— x AT * B c 

TT S 



\| 2 - 5 * k „* T ' 



x P ' 

W 00 
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FUNCTION X2D0TI T f XI ,X2 ,A) 

X2D0T=( 1.0- A) /X1*X2*X2-T*X1**(1.0-A)*X2 

°ETURN 

END 
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